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ABSTRACT 

The  correction  to  transonic  small-disturbance  theory  induced  by  the 
addition  of  contaminant  particles  (e.g.f  dust)  to  the  gas  is  analyzed.  Two- 
phase  flow  equations  governing  the  particles  and  inviscid  fluid  are  used.  The 
dusty  gas  version,  a  particular  limit  of  the  flow  equations  in  which  the 
volumetric  concentration  of  the  particles  is  small,  but  the  mass 
concentrations  of  particles  and  fluid  are  comparable,  is  employed.  In  this 
setting,  a  model  in  which  the  gas  and  particles  may  be  viewed  as  an  equivalent 
gas  with  modified  properties,  is  derived.  This  "generalized  gas"  model 
behaves  like  a  normal  gas  with  a  modified  value  of  y  (the  ratio  of  specific 
heats).  Using  this  model,  a  simple  method  of  analyzing  transonic  small- 
disturbance  theory,  by  employing  a  modified  transonic  sim:'  rarity  parameter,  is 
used  to  account  for  the  effect  of  the  dust. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  fluid  flow  problems  of  interest  concern  the  behavior  of  a  gas  which 
has  been  contaminated  with  small  particles  of  dust.  The  presence  of  the  dust 
can  cause  significant  changes  in  the  flow,  and  it  is  important  to  analyse  an 
explanation  for  this  phenomenon.  This  is  done  by  examining  a  model  in  which 
the  gas  and  dust  exchange  heat  and  momentum.  In  the  limit  of  low  vol metric 
concentrations  of  dust,  but  with  strong  coupling  between  the  phases,  the  model 
equations  are  closely  approximated  by  the  equations  for  an  adiabatic  ideal 
gas,  with  modified  values  of  the  density  and  ratio  of  specific  heats.  By  the 
use  of  similarity  transformations  of  these  equations,  it  is  possible  to  relate 
solutions  of  flow  problems  for  a  gas  with  dust  to  solutions  of  corresponding 
problems  for  a  clear  gas,  thus  giving  an  explicit  way  of  calculating  the 
effect  of  the  dust  on  the  flow.  Because  of  their  Bimple  form,  the  equations 
of  transonic  flow  are  used  to  provide  an  example  of  this  procedure.  It  is 
found  that  the  transonic  flow  around  a  thin  airfoil  for  a  gas  with  dust  is 
equivalent  to  the  flow  around  an  airfoil  with  modified  thickness,  at  a 
different  free-stream  mach  number. 
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TRANSONIC  SMALL-DISTURBANCE  THEORY  FOR  DUSTY  GASES 
Donald  A.  Drew*  and  Fredrick  J.  Zeigler** 

1.  INTRODUCTION 

The  design  of  machinery  which  utilises  or  involves  the  flow  of  a  gas  is  important  in 
many  practical  situations.  Often  such  machines  must  be  operated  in  an  environment  where 
the  gas  is  contaminated  with  small  particles  (dust)  which  may  degrade  performance,  cause 

wear  or  necessitate  filtering,  all  with  undesirable  economic  conseguences. 

In  order  to  better  understand  the  flow  of  both  the  dust  and  gas,  we  examine  a  model  in 

which  the  gas  and  dust  exchange  heat  and  momentum.  The  equations  of  conservation  of  mass, 

momentum  and  energy  for  each  material  are  simplified  by  assuming  low  volumetric 

concentrations  of  a  relatively  heavy  dust.  This  set  of  dusty  gas  equations  is  then 

further  approximated  by  assuming  strong  coupling  between  the  materials.  The  resulting 

system  of  equations  is  analogous  to  the  equations  for  the  adiabatic  motion  of  an  ideal  gas, 

except  that  the  gas  constant  Y  and  the  density  p  are  modified  to  reflect  the  heat 

m 

capacity  and  the  density  of  the  dust.  He  further  show  that  this  generalized  gas  supports 
discontinuities  in  properties  (generalized  schocks)  which  consist  of  a  shock  in  the  gas, 
followed  by  a  relaxation  back  to  equilibrium  (both  thermal  and  mechanical)  of  the  dust 
particles. 

We  note  that  the  addition  of  a  moderate  amount  of  dust  (0.1%  by  volume)  can  increase 
the  Mach  nunber  of  a  flow  by  an  appreciable  amount  (~  25%).  This  could  have  serious 
consequences  in  operation  and/or  efficiency  in  many  devices. 

Finally,  we  discuss  the  transonic  small- disturbance  theory.  It  is  shown  that  the 
effect  of  the  dust  is  to  modify  the  flow  in  such  a  way  that  it  is  equivalent  to  the  flow  of 
a  clear  gas  at  a  different  Mach  number  around  an  airfoil  of  a  different  thickness.  This 
suggests  a  way  to  study  the  behavior  of  airfoils  of  different  thicknesses  (but  similar 
shapes)  may  be  to  add  dust  to  the  flow. 
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where  2.1,  2.2,  and  2.3  are  equations  of  conservation  of  mass,  momentum,  and  energy  for  the 

particle  and  fluid  phases.  Here  q  represents  the  velocity  (p  for  particles,  f  for 

fluid),  a  is  the  dust  volumetric  density  »  dust  volume/total  volume,  p  denotes  the 

phasic  density  (for  the  particles,  p  -  dust  mass/dust  voltste),  p  is  the  pressure,  and 

P 

e  is  the  internal  energy  per  unit  mass. 

In  addition,  equations  of  state  are  needed  for  both  phases,  the  dust  is  assisted  to  be 
incompressible 

p  -  const  t  (2.4) 

P 


-2- 


for  the  fluid,  the  Ideal  gaa  law 


P  -  *  Vf 


(2.5) 


holds.  The  particles  ars  assisted  to  be  calorlcally  perfect,  so  that 

e 


<P> 


e(p)T  , 

P  v  p 


where  cv  denotes  the  constant  volume  specific  heat  of  the  particles;  and,  for  the  gas. 


(f)„ 

ef  “  cv  Tf  p/pf 


(2.6) 


The  terms  on  the  right  of  (2.2)  and  (2.3)  involving  and  reflect  momentim  and 


energy  transfer  between  the  particles  and  gas.  For  bM  we  assume  the  general  form 


w  3  ,♦  ♦  ,  1 
bM  "  e  VflqP  "  af 1  *  I 


(2.7) 


where  cD  is  the  drag  coefficient,  and  a  is  the  effective  particle  radius.  If  the 


relative  velocities  between  the  two  phases  are  small. 


cD  ■  24/Re 


(2.8) 


where 


2a0 


(2.9) 


0  being  a  reference  velocity.  Then 

bM  “  2  PfVf  (2.10) 

a 

according  to  the  stokes  law. 

Expressed  in  terms  of  the  Nusselt  number  Nu,  the  heat  transfer  coefficient  H,,  is 

(2.11) 


"h  ■  f  *  •  h 


where  k  is  the  coefficient  of  heat  conductivity.  By  including  the  Prandtl  nusber 


■  iifA,  this  may  be  written 


The  dusty  gas  limit  of  the  previous  two-phase  flow  equations  corresponds  to  the  limit 


a  <<  1,  p  <<  p  and  ap  ~  p  .  in  order  to  exhibit  this  limit  explicitly  we  first 
f  p  p  f 

scale  the  variables  and  put  the  equations  into  nondimens ional  form. 

Thus ,  we  asstmie  the  problem  contains  a  typical  length  scale  L(  velocity  scale  U, 
gas  density  scale  r ,  temperature  scale  Tg»  and  volumetric  concentration  scale  K.  He 
then  set 

x  -  x/L  (3.1a) 


t  -  tU/L 


o(x,t)  -  Aa(x,t) 


Pf(x,t)  -  Tpf(x,t) 


*  <x,t)  -  (Jq  (x,t) 
P 


V 


♦  ♦  ♦  .♦  *. 
af(x,t)  -  tJqf(x,t) 


p(x,t)  -  ru2p(x,t) 


♦  2“  ♦  * 
c  (x,t)  -  0  e  (x,t) 
P  P 


c  (x,t)  •  02e  (x,t) 
r  r  • 


T  (x,t)  -  T  T  (x,t> 
P  0  P 


Tf(x,t)  -  TQTf(x,t) 


(3.1b> 

(3.2a) 

(3.2b) 

(3.2c) 

(3. 2d) 

(3.2e) 

( 3.2f ) 

(3.2g) 

(3.2h) 

(3.21) 


The  following  dimensionless  combinations  will  be  usedt 


s  -  r/P 


e  -  ru/b^L 


h  -  HMTQ/bMU 


(3.3a) 

(3.3b) 


(3.3c) 


Y  -  R'Tjj/O 


c  -  c<p) tV 
v  0 


(3. 3e) 


;  -  c(f)Tv 

v  0 


(3.3f ) 


The  dimensionless  two-phase  flow  equations  are 


+  V«  (aq  >  -  0 


(3.4a) 


3(1  -  Aa )p 


+  7*  (1  -  Aa)pfqf  =  0 


(3.4b) 


3<to  *  ~  ; 

^  *  v7q*  *  "sVp  + «  <qf "  qP) 


(3.5a) 


(1  -  Aa)pf(-T^  +  qf-7qf)  -  -{1  -  Aa)7p  ♦  £  a(5p  -  qf) 
3 1 


.  3(c  +;qj  ?  *  *  -  a,  , 

<*[ - 2 — .  "E"  +  q  *V(e  +  2  qD ^  “ 

3t  P  P  P 

A  A  A  A  A 

A  A  A  AA.  _  *  i  ^  ^  Vs  0  A  A  A 

-  -sVa-pqp  -  sp7a*q  +  —  a(qf  -  qp)*qp  ♦  —  a(Tf  -  Tp) 


*  i  *2 

-  -  3(e,  +  r  q.)  *  .  .  ,  , 

-  Aa)pf( - : - +  V7<ef  +  2 


(3.5b) 


(3.6a) 


-7(1  -  to)-Mf  -  Ap7a»qp  -  J  a(qf  -  qp)-qp  +  f  a(Tp  -  Tf)  . 


(3.6b) 


The  dusty  gas  limit  corresponds  to  A  ♦  0,  s  ♦  0  (since  p^  <<  pp  =  const  implies 
-  s  <<  1),  with  f  »  A/S  remaining  fixed.  The  dusty  gas  formulation  is 


—  +  7« (no  )  -  0 
3t  P 


(3.7a) 


3p.  .  .  ; 

—  4-  7«rt  rt  s  0 


(3.7b) 


4.  GENERALIZED  GAS  MODEL 

Let  us  see  the  effect  which  the  effective  particle  radius,  a,  has  on  the  dusty  gas 
equations.  By  using  the  Stokes  drag  fornula,  we  may  write 


e 


2  Ua _  J_ 

9  y  AL  * 


(4.1) 


Thus,  for  sufficiently  small  particles,  e  will  be  small.  In  fact,  for  many  flows 
Re  a  and  a/L  are  small,  but  A  is  also  small.  The  rest  of  our  discussion  will  be 
limited  to  the  regime  e  <<  1. 

♦  ♦ 

For  this  case,  inspection  of  equation  3.S  shows  that  q  »  q  (we  now  drop  the  carets 

f  P 

for  convenience)  unless  the  accelerations  are  large,  ttiis  suggest  a  model  which  we  term 

♦  +• 

the  generalized  gas  model,  which  has  the  property  that  q^  *  +  0(e),  except  in  places 

where  the  flow  fields  change  rapidly.  We  shall  call  these  regions  of  rapid  change 
generalized  shocks.  As  we  shall  see,  this  model  is  analogous  to  a  gas  with  changed 
properties. 

♦  ♦  ♦ 

Consider  the  flow  fields  a,  q  ,  q  ,  etc.  as  functions  of  x,  t,  and  e.  In  the 

P  P 

generalized  gas  model,  we  consider  a  limit  of  the  dusty  gas  equations  in  which  e  ♦  0, 

♦ 

with  x,  t  held  fixed.  Away  from  generalized  shocks,  we  therefore  consider  an  expansion 

a (x,t>c )  »  a(01(x,t)  ♦  ea(1,(x,t)  *  ...  (4.2) 

with  similar  expressions  for  the  other  flow  quantities. 

Substituting  these  expressions  into  the  dusty  gas  equations  and  equating  terms  of 


equal  order  to  zero  gives  the  following: 


pm?  + vo>  ♦ 

*[— ■  *  ;*V,CT  ♦  i  «*>]  -  <%"  -  K"'-'  +  *"1"  -  Tpn) 


(0)r3(6f  +  2  q  ’  ♦  (0)  1  2.i 

P,  t - zZ - +  «,7(ef  +lq)} 


„  (0)+  (0)  +(1>  ♦  <1 ) ,  -►  .  _ _ (0  ) ,_( 1 )  _(1). 

-Vp  q  -  a  (qf  -  q^  )*q  +  ha  (T^  -  Tf  ) 


(0)  (0), 

p  *  rpf  T 


„(0)  _  *  ^  <0)  .  <0) 

ef  “  cT  +  p  /pf 


(4.11] 


(4.12) 


This  model  allows  us  to  derive  the  equations  needed  at  the  lowest  order.  Adding 


times  (4.7)  to  (4.8)  yields 


»J°’  ♦  aloln(ff  *  ;-'5) 


(4.13J 


A  similar  combination  of  (4.9)  and  (4.10)  gives 


*  ^’7cTl  +  “(0)f[f^  +  9,VcTl  + 

.  (eJ">  .  .<°>r,[!ill  . 


(0)  ,.(0) 


3  ( p  /p 


-  +  q«?(p  /p  -  -V*p  q  . 


(4. 14) 


Additionally,  (4.5)  and  (4.6)  may  be  combined  as 


3(p!0)  f0,  3p 

- vZ -  +  VMpi  +  a  f )q  -  0  »  rr^  +  7»p_q 


(4.15) 


(0)  <0) 

where  p  »  p  ♦  a  f  is  the  mixture  density.  Using  (4.5)  and  (4.6)  we  obtain 
m  r 


(0 ) 

3  f  a  \  +  0  (0) 

it  tn(loT^  +  <5*7tn(a  /pf  )  =  o  . 

pf 


(4. 16) 


0  (0 ) 

Therefore  a  /pf  is  constant  for  a  fluid  particle. 

The  energy  equation  may  be  rewritten  in  order  to  yield  further  information.  First, 
taking  the  dot  product  of  q  with  (4.13)  shows  that 


(p'0)  +  a(0)f)[— —  +  q*7(-j  q2)]  *  -q*7p 


(0) 


(4. 17) 


Then,  from  (4.14),  we  are  left  with 


(0)  (0) 

(p^>;  4.  a(o)f=)[|j 4.  5-t]  4.  p(°)[!iL-fr-l  +  ^7(pto>/P;o>,] 


«0)_  ♦ 

-p  V»q 


(4. 18) 


With  the  expressions  (4.6)  and  (4.11),  this  yields 


(Pj0)c  +  a(0,fc  +  Pf°ir>[|^  +  q*7T]  -  rTt~ 3t —  +  q*7Pf0>] 


(0) ,r3T  .  * 


3p 


(0) 


(4.19) 


Therefore,  following  a  fluid  particle. 


(0)  -f-1 

T  -  Const* (p^  V 


where 


(4.20) 


t  * 


(0), 

a  f  c 


(0)  *  _ 

Pf  c  +  r 


*  .  ^  (0).  (0)  (0), 

c  +  r  +  a  fc/pf  1  +  a _ f  _ c 


(0) 

Pf  c  +  r 


(4.21) 


is  constant  on  each  streamline. 


0  . 
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is  constant  following  a  fluid  particle 


where  PQ  and  pQ  are  constants.  Thus,  a  dusty  gas  with  snail  dust  particles  behaves 
like  a  gaa  with  a  modified  y.  We  emphasize  that  this  derivation  assumes  that  we  are  not 
near  a  shock. 

We  note  that  the  speed  of  sound  for  a  generalized  gas  is  given  by 

;  <0>  *  (0) 

2  „  d£_  la  _  :p_f_  , .  Ipf-  ;2 
dp  p_  TP  ^  (0)'  TP  ' 

HI  In  In  n 

where  a  is  the  speed  of  sound  in  the  clear  gas.  Moreover,  since  p  «  pi°^  ♦  o/°*f, 

*  (0)  .  “  f 


5.  GENERALIZED  SHOCK  STRUCTURE 

In  ord«r  to  derive  expreasions  which  measure  the  correction  to  transonic  small- 
disturbance  theory  caused  by  the  addition  of  dust  to  the  gas,  we  must  examine  how  shock 
waves  fit  into  this  framework.  Hie  generalized  gas  equations  (4.23)-(4.25)  are  not  valid 
in  regions  of  rapid  change.  In  this  section  we  shall  derive  equations  valid  for  such 
regions,  for  the  purpose  of  obtaining  shock  relations  which  can  be  used  with  the  equations 
of  the  previous  section.  We  assizne  that  the  conditions  for  a  dusty  gas  are  not  violated  in 
the  shock,  so  we  shall  start  the  analysis  by  using  (3.9)-(3.12). 

we  consider  only  planar  shocks,  with  the  generalized  shock  location  given  by 

x  *  x  (t)  ,  (5.1) 

s 

y,  z  arbitrary. 

A  coordinate  system  fixed  in  the  shock  will  be  used,  with  the  positive  x-axis  in  the 
direction  of  net  flow  through  the  shock.  Setting  X  =  x  -  xg(t),  Y  *  y,  Z  *>  z,  T»t 
transforms  the  dusty  gas  equations  to 


3a 

3T 


+ 


lit 


3  (au  )  3  (aV  )  3(aW  ) 

- B_  *  - B_  .  - E_ 


3X  + 

3Y  + 

_ 

3z 

3p  U 

3p  V 

3p  W 

*  *  f  * 

—LL  + 

f  f 

(5.2) 


(5.9) 


3w  3w.  9w  3w 

P  f -  +  U  —  4  V  4.  u  — 1  a  E  4  2  (Jj  -  W  ) 

Pf'-9T  °f  3x  f  3y  wf  3z  >  3z  e  p  V 


3t 


3t 


3t 


3T 


cf  — ^  +  n  — E  4  v  — P  4  y  — ) 

13t  p  3X  p  3Y  f  p  3Z 


!)  -  7T  tTf  •  Tp> 


(5.10) 


.  3T 


3t 


3T 


3T. 


30.  3V.  3W. 


Pfc(3T  +  Df  3x  +  Vf  3Y  +  Wf  3Z  ^  *P^3X  +  3Y  +  3Z  ^ 


-  7  ((0.  -  U  )2  +  (V,  -  V  )2  +  <w  -  W_) 2]  +  7  a(T  -  T_) 
e  f  p  f  p  fp  epf 


(5.11) 


dx 

In  th«  above,  O  *  u  .  -  — V  ,  -  v  ,  W-  w  .  The  choice  of 
p,f  p,f  dt  P»r  P»f  P#“  P»f 

coordinate  system  dictates  U  „  >  0. 

P,f 

To  observe  the  details  of  the  shock  structure,  we  examine  a  distinguished  limit  of 
these  equations  most  appropriate  for  regions  of  rapid  change.  In  order  to  get  the  richest 
equations,  we  choose  £  «  X/e .  Then  we  have 


3  (au  ) 


3(Pf0f) 

H 

30 


"pH  *?"f-  V 


3V 

up  -  7  <vf  “  V 


3W 

UP  3£*  “  f  <Mf  •  V 

»u  9 

PfWf  3£  ‘  H  +  0(UP  •  Uf> 


3V 

P-0.  TT*  -  a  (V  -  V,) 
f  f  3£  p  f 


3H 

P,U  TJ*-  -  U(M  -  »  ) 

f  f  3£  p  f 
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(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 

(5.19) 


*4 


-  ■  -1 
.  "J 


. '  J 

.  i 


5 


p  at  i  t  p 


*  3Tf 

pfcDf  IT 


3a. 


-p  ■TT’  -  «t<0,  -  0  )2  ♦  (V.  ♦  V  )2  ♦  (W_  -  V*  )2]  +  ha(T  -  T .) 
of.  Ip  Ip  f  p  p  I 


Boundary  conditions  are  given  as: 


dx 


U.,0  ♦  u"  -  —)  v  ,v  ♦  V ±J  w  ,w  ♦  *t, 
f  p  at  r  p  f  p 

±  ±  + 

T-,T  ♦  T»  a  ♦  a  »  P,  *Pl  as  5  ♦  • 

f  p  r  r 


(5.21) 


Let  us  consider  the  possibility  that  these  eauations  have  a  shock  solution,  that  is,  a 
discontinuous  solution  which  satisfies  appropriate  jump  conditions.  Equations  (5.12)- 
(5.20)  give  rise  to  the  jump  conditions 


too  ]  -  0 
p 

(5.22) 

[pfuf3  -  o 

(5.23) 

[U2]  -  0 
p 

(5.24) 

[U  V  ]  -  0 

(5.25) 

p  p 

[0  W  ]  -  0 

(5.26) 

p  p 

lpf°f  +  p]  ■  0 

(5.27) 

tpfu  V  -  o 

(5.28) 

(PfUfWfl  -  0 

(5.29) 

[COT  ]  -0 

(5.30) 

p  p 

tPf<=UfTf  +  p0f]  -  0 

(5.31) 

! 


-13- 


Condition*  (5.22),  (5.24),  (5.25),  (5. 26),  and  (5.30)  show  that  tha  particl* 
proparti**  ar*  continuous  across  a  shock.  Tha  other  condition*  show  that  tha  ga* 


(5.36) 


To  get  the  total  momentum,  we  add  at  times  (5.14)  to  (5.17),  which  yields 

2  2 

afU  *  p.U.  ♦  p  ■  g(t)  . 

P  11 

By  an  argument  similar  to  the  one  leading  to  (5.35)  we  find 


p«tu" 


dx  2 


dx  2 


+  p 


(5.37) 


By  talcing  a  linear  combination  of  (5.20)  and  (5.21),  and  using  (5.14)  and  (5.17)  on 

the  resulting  expression,  we  may  derive  the  energy  equation 

*  3  3 

OfcUpTp  ♦  PfcUfTf  ♦  OfUp/2  ♦  Pfuf/2  +  P°f  "  e(t>  •  (5.38) 

The  corresponding  jump  condition  is 


dx 

dx  3 

dx 

pfc>( 

>■  - 

+  p"(u"  •  5*) 

dx 

dx 

(o+fc 

+*  /  + 

+  pfc)(u  - 

.  1  ♦/  +  _ S\  3  .  +f  + 

+  2P«lU  ‘  dt  )  +p(U 

dx 


(5.39) 


Two  equations  cannot  be  explicitly  integrated  across  the  generalised  shock.  They  can 
be  taken  to  be 

*.  i 

(5.40) 


Up  "  f  tUf  -  °p>  ' 


dT 


CU  -  “  ( T  -  T  )  . 

p  dt  f  f  p 


(5.41) 


we  note  that  (5.40)  and  (5.41)  govern  the  non-equilibrium  evolution  of  the  velocities 

and  temperatures  hack  to  the  equilibrium  values  U+  and  T+.  For  the  generalised  shock 

structure  we  describe,  it  is  sufficient  that  these  equations  ing>ly  U  ♦  U  and  T  ♦  T 

p  f  p  f 

as  5  ♦ 

The  junp  conditions  (5.35),  (5.37),  and  (5.39)  show  that  for  a  generalised  gas,  the 
generalized  shock  relations  are  precisely  those  for  a  gas  with  equivalent  properties.  The 


shock  thickness  if  0(e  ) 


6.  SMALL  DISTRPBANCE  THEORY 

With  the  derivation  of  the  generalized  gas  model  (4.23)-(4.25)  and  generalized  shock 
conditions  (5.35),  (5.37),  and  (5.39),  we  have  achieved  considerable  simplification  of  the 
original  problem.  The  dusty  gas  may  be  regarded  as  an  equivalent  gas  with  a  modified  value 
of  Y>  Ass  lining  that  we  have  potential  flow,  it  is  therefore  possible  to  formulate  a 
small-disturbance  theory  for  transonic  thin  airfoils  in  dusty  gases,  analogous  to  that  in 
[2] .  The  transonic  similarity  parameter  for  this  case  is  a  function  of  the  new  value  of 
Y.  It  will  be  shown  how  this  modified  similarity  parameter  is  related  to  the  usual  one 
for  the  case  of  no  dust,  so  that  a  method  of  estimating  the  effect  of  the  dust  on  the  usual 
small-disturbance  theory  may  be  achieved. 

We  formulate  the  boundary  value  problem  for  a  thin  airfoil,  in  a  dusty  gas  and 
travelling  in  the  transonic  range,  as  follows 


Figure  6-1 

As  shown  in  Figure  6-1,  the  free-stream  velocity  is  0  (in  the  x  direction  only), 
where  the  coordinates  x,  y  have  been  normalized  with  respect  to  the  airfoil  chord.  The 
function  F  which  defines  the  airfoil  surface  satisfies  max  |F  (x)  -  F  ( x) I  -  1. 

xeto.i]  u  1 

we  assisne  that  a  potential  9  exists,  so  that  the  flow  components  are  given  by 

79  -  (U  ♦  u,v)  -  Ue  ♦  q  (6.1) 

♦ 

where  u,v,  represent  the  distrubance  velocities,  and  e^  is  the  unit  vector  in  the  x 


direction 


The  potential  4  satisfies  the  system 


>♦ 


*2  2 

2H  ♦  ♦  (a  -  •  )♦ 

x  y  xy  y  yy 


(6.2) 


(6.3) 


Hare,  a  is  the  local  speed  of  sound  in  the  dusty  gas,  which  we  are  going  to  distinguish 
fro*  a,  the  sound  speed  in  the  gas  without  dust. 

6  satisfies  the  boundary  condition 


4  (x,4f  .<x)> 

-1 - HJ: -  .  Sr.  (x) 

4  (x,6p  ,<x)>  0 


(6.4) 


X  •  u,t' 

of  tangent  flow  along  the  wing,  and 

4  ♦  U  as  x  ♦  —  ,  (6.5) 

x 

corresponding  to  unifora  flow  for  upstream.  He  also  assume  that  the  Kutta  condition 
holds:  the  flow  leaves  the  trailing  edge  smoothly. 

Under  the  condition  of  small  distrubances ,  we  may  derive  a  simpler,  approximate 
equation  for  the  problem,  this  derivation  is  achieved  through  a  limit  process  expansion 

~  1/3  •  1  “ 

for  ♦,  based  on  the  limit  process  4*0,  with  x,y  -  6  y,  and  K  -  — irTT-  fixed  as 


,2/3 


M  *  1.  It  has  the  form 


then  4  satisfies 


4(x,y»6)  -  U{x  ♦  62/34(x,y)  *  ...} 


(K  -  (y  ♦  1)4  >4  ♦  4—  -  0 

x  xx  yy 


and  the  boundary  conditions 


4~(x,0  )  «  p*  <x)  0  <  x  <  1 

y  Ugfc 

4  ,4~  *0  as  x  »  -m  . 
x  y 

In  this  context, the  Kutta  condition  may  be  written  as 

1VtE  ■  °  • 


(6.6) 

(6.7) 

(6.8) 

(6.9) 

(6.10) 


where  TE  means  trailing  edge. 


Additionally,  shock  juap  conditions  must  bs  imposed  in  ordsr  to  ha vs  a  couplets 
problea  for  ♦  .  the  conditions  may  bs  darivsd  f roa  ths  continuity  of  ths  transonic  asss 
flux  vector 

Pq  -  P.0(sx(1  +  «4/3(K*x  -  *-*-  *2>  ♦  ...)  ♦  +  ...  .  (6.11) 

If  ths  shock  surface  is  given  by 

S(x,y)  -  x  -  g(y)  -  0  , 


ths  jtaap  condition  is 

Wx'^-pOs*  "  0  *  <«•«> 

Furthermore,  there  is  no  juap  in  the  tangential  component  of  velocity  across  the  shock, 
which  may  be  guaranteed  by  imposing 


-  0  .  (6.13) 

The  fora  of  (6.7)  shows  that  for  two  flows  at  different  transonic  Mach  nissbers,  but 

with  the  same  values  of  K,  these  two  flows  will  be  geoaetrically  similar  (the  difference 

will  be  that  the  'size'  of  the  disturbance  will  be  determined  by  different  factors  of 
2/3 

4  ) .  An  analogous  similarity  law  holds  for  the  gas  alone,  governed  by  the  similarity 

paraaeter  K.  Both  of  these  similarity  laws  relate  families  of  flows,  for  different  values 
of  the  displacement  thickness  at  a  corresponding  Mach  number. 

We  wish  to  see  how  the  change  in  the  ratio  of  the  specific  heats  for  the  dusty  gas, 

A 

Y,  affects  the  flow.  This  will  be  done  by  finding  a  correspondence  between  the 
similarity  parameters  in  the  two  cases. 

We  shall  develop  the  correspondence  by  transforming  the  boundary  value  problea  (6.2)- 
(6.9)  for  the  dusty  gas  into  a  problem  for  the  gas  without  dust. 

Let  ♦  (*)  ■  ♦  be  the  solution  operator  for  the  dusty  gas  boundary  value  problea.  By 
defining 

i  -  w"2K  (6.14a) 

y  »  «y  “  j^y  (6.14b) 

♦  -  u>*  (6.14c) 


where  w 


f*  ±  llV3 

'Y  *  1  ' 


and  6 


3. 

u  5 , 


we  see  that  (6.7),  (6.8)  become 
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(6.15) 


(X  -  <Y  ♦  D*  )♦  +  -  0  , 

x  xx  Tyy 


Vx,°  ’  "  W"’  <><*<'• 


(6.16) 

The  other  boundary  conditions,  (6.9)  and  (6.10),  transform  readily  in  the  correct  manner. 
The  shock  jump  condition  (6.12),  since  it  can  be  derived  from  the  divergence  form  of  (6.7), 
iipoeea  no  additional  restriction  upon  transformation  by  (6.14).  This  is  also  true  for 
(6.13). 

aw  * 

Therefore,  under  (6.14)  the  entire  problem  for  +(x,y>K)  transforms  to  the 
corresponding  problem  for  £(x,ysic),  which  is  the  case  of  the  gas  without  dust.  Comparing 
(6.7)  and  (6.15),  we  conclude 

*(K)  -  .  (6.17) 

Now,  the  similarity  parameter  X  for  the  new  problem  is  of  the  form 


i  -5! 


i-5! 


i-5! 


K  - 


y2/3 


,  3,, 2/3  2,2/3 

(w  8)  u  8 


-2* 
U  K 


(6.18) 


from  which  we  conclude  Combining  (6.17)  and  (6.18)  gives 


♦  (X)  -  u<i(u"2k  )  . 


(6.19) 


This  relation  is  equivalent  to 


~  *  -  -2- 
♦  (x,y»x)  -  u4(x,y»w  X)  . 


(6.20) 


To  complete  the  details,  we  note 

.  i-m! 


,2/3 


where  *  « 


t(i 


.  CO  > 


from  Section  4,  so  that 


t(i 


(0). 


.  <0» 


1  - 


,2/3 


(6.21) 


where  y  is  given  by  (4.21). 
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Thua,  we  see  from  (6.20)  that  the  small  disturbance  problem  for  a  dusty  gas  is 
equivalent  to  one  for  a  clear  gas,  with  a  modified  amplitude  of  disturbances,  a  changed 
wing  thickness  ratio  and  similarity  parameter. 

(0)  (0) 

The  transformation  parameter  u,  as  a  function  of  a  f/pf  ,  is  shown  in  Figure 
6.2  for  c/(c  +  r)  “  1.0,  which  is  a  representative  value  for  many  materials. 


Figure  6.2.  w  vs.  a  f/pf 

Thus,  for  example,  a  dusty  flow  with  a  mass  loading  of  2.0  corresponds  to  a  value  of 
w  of  1.21,  which  means  that  if  the  airfoil  has  thickness  6,  the  flow  is  equivalent  to 
that  of  a  clear  gas  around  an  airfoil  of  the  same  shape,  but  of  thickness 

-  3  _ 

o  «  w  6  »  1.776.  The  equivalent  similarity  parameter  K  is  K  ■  g  l  •  0.68K.  This 
implies  that  for  a  given  airfoil,  of  thickness  60»  the  solution  in  a  clear  gas  for 
thickness  6t  -  u360  at  similarity  parameter  K,  can  be  obtained  from  the  solution  in  a 
dusty  gas  with  a  value  of  id  equal  to  u  -  (61/6(J)1/3,  and  similarity  parameter 
■  oi2K,.  TOe  loading 


can  be  obtained  from 


■o  that 


<W 


1/3 


+  1x1/3 
+  V 

+  1) 


Y 
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